High-order harmonic generation enhanced by XUV hght 
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The combination of high-order harmonic generation (HHG) with resonant xuv excitation of a core electron into 
the transient valence vacancy that is created in the course of the HHG process is investigated theoretically. In 
this setup, the first electron performs a HHG three-step process whereas, the second electron Rabi flops between 
the core and the valence vacancy. The modified HHG spectrum due to recombination with the valence and the 
core is determined and analyzed for krypton on the 3(i — >■ 4p resonance in the ion. We assume an 800 nm laser 
with an intensity of about 10^* and xuv radiation from the Free Electron Laser in Hamburg (FLASH) 

with an intensity in the range 10^^-10^^ prediction opens perspectives for nonlinear xuv physics, 

attosecond x rays, and HHG-bascd spectroscopy involving core orbitals. © 2011 Optical Society of America 
OCIS codes: 190.2620, 140.2600, 190.7220, 260.6048. 



High-order harmonic generation (HHG) by atoms in 
intense optical laser fields is a fascinating phenomenon 
and a versatile tool; it has spawned the field of atto- 
science, is used for spectroscopy, and serves as a light 
source in many optical laboratories [T] . Present-day the- 
ory of HHG largely gravitates around the single-active 
electron (SAE) approximation and the restriction to 
HHG from valence electrons [iHl]- 

Several extensions to the SAE view of HHG have been 
investigated previously. A two-electron scheme was con- 
sidered that uses sequential double ionization by an op- 
tical laser with a subsequent nonsequential double re- 
combination; in helium this leads to a second plateau 
with about 12 orders of magnitude lower yield than 
the primary HHG plateau [5]. Two-color HHG (opti- 
cal plus xuv light) has been studied in a one-electron 
model [6] and with many-electron effects included by a 
frequency-dependent polarization [7]; the xuv radiation 
assists thereby in the ionization process leading to an 
overall increased yield [B] and the emergence of a new 
plateau [7], the latter, however, at a much lower yield. 
The above schemes suffer from tiny conversion efficiency 
beyond the conventional HHG cutoff (maximum photon 
energy). 

We propose an efficient two-electron scheme for a 
HHG process manipulated by intense XUV light from the 
newly constructed free electron lasers (EEL) — e.g., the 
Free Electron Laser in Hamburg (EL ASH). Our prin- 
cipal idea is sketched in Eig. [1] In the parlance of the 
three-step model IH[3], HHG proceeds as follows: (a) the 
atomic valence is tunnel ionized; (b) the liberated elec- 
tron propagates freely in the electric field of the optical 
laser; (c) the direction of the optical laser field is reversed 
and the electron is driven back to the ion and eventu- 
ally recombines with it emitting HHG radiation. The 
excursion time of the electron from the ion is approxi- 
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Eig. 1. (Color online) Schematic of the three-step model 
for the HHG process augmented by xuv excitation of a 
core electron. 



mately 1 fs for typical 800 nm optical laser light. During 
this time, one can manipulate the ion such that the re- 
turning electron sees the altered ion as depicted in Eig.[T] 
Then, the emitted HHG radiation bears the signature of 
the change. Perfectly suited for this modification dur- 
ing the propagation step is xuv excitation of an inner- 
shell electron into the valence shell. The recombination 
of the returning electron with the core hole leads to a 
large increase of the energy of the emitted HHG light as 
the energy of the xuv photons wx is added shifting the 
HHG spectrum towards higher energies. A prerequisite 
for this to work certainly is that the core hole is not too 
short lived, i.e., it should not decay before the continuum 
electron returns. 

The spatial one-electron states of relevance to the 
problem are the valence state | a ) and the core state | c ) 
of the closed-shell atom. In strong-field approxima- 
tion, continuum electrons are described by free-electron 
states I fc ) for all fc S M"^ U . The associated level en- 

ergies arc Ea, Ec, and We need to consider three 
different classes of two-electron basis states to describe 
the two-electron dynamics: first, the ground state of 
the two-electron system is given by the Hartree prod- 
uct I a ) dg) I c ) ; second, the valence- ionized states with one 
electron in the continuum and one electron in the core 
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state are | fc ) (g) | c ) ; third, the core- ionized states with one 
electron in the continuum and one electron in the valence 
state are \k) \a). We apply the three assumptions of 
Lewenstein et al. [4] in a somewhat modified way by 
considering also phenomenological decay widths of the 
above three state: Tq and Fq to account for losses due to 
ionization by the optical and xuv light for | a ) Cg) | c ) and 
\k) ig) |c), respectively, and Tc to represent losses from 
ionization by the optical and xuv light and Auger decay 
of core holes for | fc ) | a ) with radiative decay of the 
core hole being safely neglected. Further, the xuv light 
induces Rabi flopping in the two- level system of | fc ) ® | a ) 
and \k) | c). 

The two-electron Hamiltonian of the atom in two-color 
light (optical laser and xuv) reads H ~ Ha + Hh + -f^x; 
it consist of three parts: the atomic electronic struc- 
ture Ha, the interaction with the optical laser Hi^, and 
the interaction with the xuv light Hx- We construct H 
mostly from tcnsorial products of the corresponding one- 
particle Hamiltonians Ha, hi^, and hx- The interaction 
with the optical and xuv light is treated in dipole ap- 
proximation in length form [8|. 

We make the following ansatz for the two-electron 
wavepacket (in atomic units) 

\^,t) = a(i)c-i(^"+^'=-"^'*+'^f'*|a)(8|c) (1 

-I- / [6a(fc,t)e-5(£;a + Sc-^x)t+iIpt ^ |c) 

-t-6c(fc,t)e-^(So+B.+^x)t+iipt dS/c 

where we introduce a global phase factor based on Ip = 
—Ea + |. The detuning of the xuv photon energy from 
the energy difference of the two ionic levels \b 5 = Ea — 
Ec—'jJx- The index on the amplitudes ha{k, t) and bc{k, t) 
indicates which orbital contains the hole. 

We insert \ ,t) into the time-dependent Schrodinger 
equation and project onto the three classes of basis states 
which yields equations of motion (EOMs) for the in- 
volved coefficients. We obtain the following EOM for the 
ground-state population 



We change to the basis of xuv-dressed states using the 
eigenvectors U and eigenvalues A+, A_ of R. 

To determine the HHG spectrum, we solve Eq. ^ 
by neglecting the second term on the right-hand side 
as in Ref. — its influence is included in Fq, Fq. and 
Fc — and a constant xuv flux starting at t = and 
ending at t = Tp. The HHG spectrum is given by the 
Fourier transform of ( 'fo, ^ | ^ | i ), where D is the 
two-electron electric dipole operator [8], |^'o,i) is the 
ground-state part of the wavepacket ([T|) and |vl'c7^) is 
the continuum part, i.e.. 
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Here, w = U ^ (q) , the ponderomotive potential of the 
optical laser is C/p, and Foj(t), C(r) are defined as in 
Ref. [4] augmented by cjl and with /p replaced by ^ -|- 
Xp. Further, Jjv arc Bcssel functions and bM-N,i{T) are 
the coefficients defined in Eq. (17) of Ref. [4] for valence- 
)and core-hole recombination. Further, 
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Neglecting the factor e 2 for now, we see that 
the hM,N.i{^,T) peak at = ^M,N,i ~ (2 (il/ + N) + 
<5i,a) wl + Si^ctdx- In other words, for i = a the peaks are 
at the positions of the harmonics from optical laser-only 
HHG; for i ~ c, the harmonics are shifted by lux with 
respect to the harmonics for i — a such that, in general, 
none of them coincides with harmonics from i = a. The 
harmonic photon number spectrum (HPNS) for a single 
atom — the probability to find a photon with specified 
energy — along the x axis is given by 
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ba{k,t) {a\hi^\k) d^fc. (2) 
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The other two EOMs are written as a vector equation, 
defining the amplitudes b{k,t) — {ba{k,t),bc{k,t))^ , 
the Rabi frequency i?ox for continuous-wave (CW) 
xuv light and the Rabi matrix 
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This yields for a CW optical laser electric field E\^{t) 
oscillating with angular frequency wl: 
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dt 



b{k,t) 



(R+ik^ + 2Ip)l.)b{k,t) + Ei^{t) (4) 
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with the density of free-photon states q{^) [5] and the 
solid angle fig. 

We apply our theory to krypton atoms. The energy 
levels are Ea = -14.0 eV for Kr4p [10] and E^ = 
— 96.6eV for KrSrf with a radial dipole transition ma- 
trix element of 0.206 Bohr [5]. The xuv fight has the 
photon energy ujx — Ea — Ec- The optical laser in- 
tensity is set to /oL = 3 X 10^"* at a wavelength 
of 800 nm. Both xuv and optical light have a pulse dura- 
tion of Tp = 3 ^ . The experimental value for the decay 
width of Kr 'id vacancies is Foxpt = 88 meV [TT] . The 
decay widths due to xuv ionization of the atom and 
the ion are obtained from the sum of the photoioniza- 
tion cross sections of the energetically accessible elec- 
trons [l2l[T3]; the ionization due to the optical laser is 
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Harmonic order h 

Fig. 2. (Color online) Photon number of hth harmonic 
order for xuv intensities of (a) lO^^ and (b) lO^^ 
The black solid lines show the contribution from recom- 
bination with a valence hole whereas the red dashed 
lines correspond to recombination with a core hole. The 
lines represent harmonic strengths obtained by integrat- 
ing over peaks in the HPNS. 



determined with Ref. |14] and is larger than the width for 
xuv ionization and Auger decay for the chosen param- 
eters. We find for an xuv-intensity of /gx = 10^^ 
To = 280 meV, = 1.5 meV, and = 88meV, and 
for an xuv-intensity of /ox = 10^^ Tq = 450 meV, 
Ta = 170 meV, and = 300 meV. 

In Fig.[2]we show the single-atom HPNS of HHG which 
is modified by xuv light. We find that the xuv excita- 
tion leads to two plateaus, one from valence- and one 
from core-hole recombination which overlap slightly. The 
width of the overlap can be tuned by changing the opti- 
cal laser intensity and interferences between both terms 
may occur. Even in the case of a moderate xuv in- 
tensity [Fig. ^^], the emission rate of HHG from core- 
hole recombination is substantial. This prediction can 
be interpreted in terms of an excitation (or even Rabi 
flopping) [9] of the remaining electron after tunnel ion- 
ization of an atom in the HHG process [Fig. [Ud] . This 
electron sits in a two-level system where the electron 
is either in the valence or the core level. The strength 
of the HHG emission due to core-hole recombination is 
roughly proportional to the population of the upper state 
around 1 fs. For 10^"^ the population is ~ 0.001 while 
it is ~ 0.6 for 10^^ Finally, we need to realize that 
the dipole matrix element for a recombination with a 
Kr 3d hole is not substantially different from the one 
with a Kr 4p hole thus explaining the similar yield of 
both contributions in Fig. [21d. 

The xuv radiation from present-day FELs is gener- 
ated with the self-amplification of spontaneous emis- 
sion (SASE) principle; it is fully transversally coherent 
but exhibits only limited longitudinal coherence. We find 
that the fluctuating phase does not destroy the spec- 
tra [8]. Thus our interpretation of Fig. [2] is not invali- 



dated when we relax the view assumed so far of entirely 
coherent xuv light with constant amplitude. 

In conclusion, we predict HHG light from resonant ex- 
citation of transient ions in HHG that allows insights 
into the physics of core electrons and has various ap- 
plications: it allows one to generate isolated attosecond 
x-ray pulses by ionizing atoms near the crests of a single- 
cycle optical laser pulse and selecting the highest photon 
energies by filtering [15]. This complements FEL-based 
strategies to generate attosecond x rays (Ref. [TB] and 
References therein). Our scheme has the advantage that 
the attosecond pulses have a defined phase-relation to 
the optical laser and that it can be employed at any FEL 
with moderate cost and minimal impact on other exper- 
iments. Further, our scheme has the potential to become 
an in situ probe of the dynamics of cations in strong op- 
tical fields interacting with intense xuv light. Namely, 
the HHG spectra depend sensitively on the xuv pulse 
shape; a reconstruction with frequency resolved optical 
gating (FROG) [T^ may be possible — thus offering the 
long-sought pulse characterization for SASE xuv light — 
but requires further theoretical investigation. Addition- 
ally, the emitted upshifted light due to core recombina- 
tion bears the signature of the core orbital; thus it can be 
used for ultrafast time-dependent chemical imaging |18j 
involving inner shells which is not feasible so far. This 
allows one to extend such HHG-based methods to all or- 
bitals that couple to the transient valance vacancy by 
suitably tuned xuv light. Our findings are not restricted 
to krypton but HHG spectra for resonant excitation of 
Is electrons in neon were successfully computed and will 
be discussed in future work. 
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